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\l Abstract. We analyse the propagation of particles in the nar- 

row funnel of a thick accretion disk. It is assumed that: (1) 

... 

~ J the funnel walls emit black body radiation with temperature 
decreasing outwards; (2) the magnetic and electric fields are 
bJOlongitudinal in the funnel. Such a scenario is consistent with 
^ models in which a large potential drop is induced by a rotat- 
ing massive black hole threaded by a magnetic field (Blandford 
fc Znajek 1977, Macdonald fc Thorne 1982). The interaction 
of relativistic protons with thermal photons from the funnel 
results in direct production of pairs and/or pions. We dis- 
cuss the relative importance of these processes for different 
r- H conditions in the funnel (temperature and electric field pro- 
files). Injected pairs interact with thermal photons in the 
Thomson or Klein-Nishina regime. Under some conditions, a 
^— 2 pair avalanche results, which we assume saturates in the pro- 
duction of stable bunches, containing almost equal numbers of 
electrons and positrons. As a consequence, highly collimated 
C 7-ray photons are produced. We obtain 7-ray spectra from our 
O |mnHpl in order to test its applicability to the 7-ray emitting 
■ AGNs recently reported by EGRET. 
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1. Introduction 

The recent discovery of high energy 7-ray emission (> 100 
MeV) from the class of active galactic nuclei (AGNs), known 
as 'blazars', has demonstrated unambiguously that high energy 
processes play a fundamental role in AGNs. Other types of 
AGN, which have previously been reported as sources of soft 7- 
rays (e.g. the Seyfert galaxies: MGC8-11-11, NGC 4151) have 
not been detected above 1 MeV by the SIGMA telescope or 
the OSSE and COMPTEL telescopes of the Compton Gamma 
Ray Observatory (Jourdain et al. 1992; Maisack et al. 1993), 
which can be interpreted as evidence for the variability of the 
high energy emission of these sources. 
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Models of 7-ray production in AGNs can be characterised 
by the geometry they assume for the central region (spher- 
ically symmetric or jet-like being the most popular choices) 
and by the basic mechanism of 7-ray production: hadronic or 
leptonic. The older hadronic models postulate the accelera- 
tion of protons (usually shock acceleration is suggested) up 
to extremely high energies in an isotropic picture. The inter- 
action of protons with matter and/or radiation then leads to 
the production of pairs (e.g., Sikora et al. 1987), neutrons 
(e.g., Kirk fc Mastichiadis 1989; Sikora et al. 1989) and pi- 
ons (e.g., Morrison et al. 1984; Bednarek fc Calvani 1991). 
Energetic secondary pairs and neutrons subsequently par- 
ticipate in the generation of 7-rays. However, in the case of 
blazars, the discovery of significant variability on a time scale 
of days, coupled with extremely high fluxes of 7-rays leads to 
severe problems with the energy budget in isotropic models, 
which require luminosities ~ 10 48 ergs -1 . Following the work 
of Marscher (1980), Konigl (1981), Reynolds (1982), Ghis- 
ellini et al (1985) and Ghisellini fc Maraschi (1985), Maraschi 
et al (1992) and Zdziarski fc Krolik (1993) have suggested 
anisotropic models in which the emission results from the syn- 
chrotron self-Compton mechanism operating in electron blobs 
moving with relativistic speed. Other similar models suggest 
the comptonisation by such a blob of soft radiation from the 
accretion disk (e.g., Melia fc Konigl 1989; Dermer et al. 1992, 
Dermer fc Schlickeiser 1993; Sikora et al. 1993). Recently, 
anisotropic hadronic models have also been proposed in which 
the 7-rays originate in the inner jet as a result of relativis- 
tic protons colliding with soft photons of synchrotron origin 
(Mannheim et al 1991, Mannheim fc Biermann 1992) or with 
matter of a thick accretion disk (Bednarek 1993). 

In this paper we examine a new anisotropic picture for the 
high energy phenomena observed in AGNs, in which the in- 
teraction and acceleration of electrons/positrons occurs deep 
in the narrow funnel of a thick accretion disk surrounding the 
central black hole. The inner parts of the funnel contain lit- 
tle matter, because the magnetic field lines there thread the 
event horizon of the black hole and effectively exclude mater- 
ial from the accretion disk. However, the region is permeated 
by radiation emitted by the walls. This we expect to be of 
approximately Planckian spectrum and to have a temperature 
near the base of the funnel significantly higher than that indi- 
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cated by the observations of the UV bump, but consistent with 
observations of soft X-ray excesses in some AGNs. The main 
difference between the scenario investigated here and previous 
work is that we assume charged particles to be accelerated by 
a component of the electric field parallel to the magnetic field. 
Such a configuration is possible only where the assumption 
used in ideal MHD of infinite conductivity breaks down, which 
can occur locally at sites of magnetic reconnection, or, globally, 
if insufficient charge is available to screen out the electric field 
induced by a rotating magnetic field. Most of our calculations 
apply equally well to each of these cases. But, where necessary, 
we concentrate in this paper on the latter picture and assume 
that both the electric and magnetic fields lie along the axis of 
the funnel, which contains no background plasma. Henri and 
Pelletier (1991) have developed a model using a similar over- 
all picture and propose the outer parts of the funnel contain 
a subrelativistic jet. However, in contrast to our model, they 
assume the inner core of the funnel contains a pair plasma of 
high enough density to short out the induced electric field. 

The acceleration of protons by the electric field in our sce- 
nario is limited either by energy losses in direct pair pro- 
duction or by losses in pion production on photons, depending 
on the parameters in the funnel. The acceleration of electrons 
or positrons, on the other hand, is limited by inverse Compton 
losses which can efficiently convert energy extracted from the 
electric field into high energy radiation. Only for very strong 
electric fields and low radiation temperatures does the process 
of triplet pair production take over as the major energy loss 
mechanism. Synchrotron radiation, which would play a role if 
particles were injected isotropically, is strongly suppressed be- 
cause pitch angles remain small (~ I/7) throughout the cas- 
cade. The details of the physical scenario and the description of 
the interaction of protons and electrons with thermal photons 
are discussed in Section 2. We find that the comptonisation of 
thermal radiation by electrons in the Klein- Nishina regime re- 
sults in the formation of bunches of pair plasma. Such bunches 
move through the funnel and scatter thermal photons into the 
high energy 7-ray range in a way reminiscent of the pulsar 
model described by Michel (1991). In terms of the bunch sce- 
nario, the 7-ray emission from two blazars the first with a hard 
and the second with a soft 7-ray spectrum (QSO 1633+382 
and CTA 102, respectively), is modelled in Section 3. Finally, 
in Section 4 we summarise and discuss our results. 

2. High energy radiation from the funnel of an accre- 
tion disk 

We start from a picture of the central region of an AGN in 
which a geometrically thick, radiation supported accretion disk 
exists around a supermassive black hole (for a review see Begel- 
man et al. 1984). Studies of the structure of such a disk pos- 
tulate the existence of a narrow funnel on the rotation axis. 
The funnel is filled with thermal radiation (we assume a black 
body distribution) which is emitted by the disk walls (e.g., 
Jaroszyriski et al. 1980, Abramowicz et al. 1980, Madau 1988). 
This picture of the central region of an AGN is consistent with 
observations of a UV bump in the spectrum of some sources 
(Shields 1978; Malkan fc Sargent 1982). In addition, there is 
observational evidence of a soft X-ray excess in some AGNs 
(e.g., Yaqoob fc Warwick 1991; Puchnarewicz et al. 1992 and 
references therein). A possible interpretation of this excess is 
that it comes from the inner part of the funnel where the char- 



acteristic temperature is significantly higher than that of the 
region which gives rise to the UV bump. In some AGNs - in 
particular in blazars - this component may be reduced be- 
cause of the narrow opening angle of the funnel; also, it may 
simply be too difficult to detect this component in AGNs at 
high redshifts, because it will be shifted to lower photon ener- 
gies and absorbed en route. On the other hand, some sources 
which show no evidence of beaming do have a soft X-ray ex- 
cess. Clearly, observations are inconsistent with the hypothesis 
of a unique opening angle for the funnels of all AGN. Never- 
theless, a significant temperature gradient is indicated in the 
disk funnel, and we incorporate this into our model. A similar 
approach has been taken by Madau (1988), who adopts a tem- 
perature profile ranging between a few times 10 6 K in the inner 
regions and ~ 10 K in outer regions of the funnel in order to fit 
the UV - soft X-ray excess seen in the spectra of some AGNs. 

Inside the funnel, we assume that particles are accelerated 
in a constant electric field. Large potential drops can be in- 
duced by a supermassive black hole threaded by a magnetic 
field predominantly parallel to its rotational axis (Blandford fc 
Znajek 1977; Macdonald fc Thorne 1982, for reviews see Begel- 
man et al. 1984; Rees 1984). This field, which is maintained by 
currents flowing in the funnel walls can be quite strong (of the 
order 10 4 G) and prevents plasma from the disk or funnel walls 
from penetrating into the inner parts of the funnel where the 
magnetic field lines are linked to the horizon of the black hole. 
As a result, this region of the funnel may possess insufficient 
charge to screen out the induced electric field. If only small 
amount of primary plasma is injected into this region of the 
funnel, (e.g., pairs from 7-7 annihilation or e-p plasma from 
the decay of low energy neutrons) the injected particles will be 
accelerated by the electric field. 

A large scale electric field can also be induced by an accre- 
tion disk rotating in a perpendicular magnetic field (see e.g., 
Lovelace 1976; Blandford 1976). However, in this case the mag- 
netic field lines pass through the inner disk, so that plasma is 
able to enter and may screen out the induced electric field. Al- 
ternatively, it is also possible that smaller scale electric fields 
are formed inside a disk or a jet, at points where magnetic 
reconnection takes place (e.g., Biskamp 1989). Although it is 
difficult to make a quantitative calculation, estimates from one 
recent model suggest electric fields of the order of 3.10 V cm~ 
(Haswell et al 1992). The discussion below of the equilibrium 
Lorentz factor achieved in an electric field applies to this pic- 
ture too, since the minimum size required in order for the pair 
production avalanche to be important is extremely small. How- 
ever, the degree of collimation of the accelerated particles in 
the reconnection picture is uncertain, and we base the calcu- 
lations in this paper on the case of a large scale field aligned 
along the funnel axis. 

We consider, then, the propagation of particles in the nar- 
row funnel of an accretion disk with longitudinal electric field 
E (measured in V cm -1 ) and longitudinal magnetic field. The 
orientation of the vectors of angular momentum of the disk 
and magnetic field in the funnel is chosen to be antiparallel, 
allowing the acceleration of particles with positive charge out- 
wards, although this is important only for proton acceleration. 
We suppose that there exists a mechanism of particle injection 
into the electric field deep in the funnel and discuss the fate of 
injected protons and leptons in turn. 
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2.1. Injection of protons 

The important interactions of protons with soft photons are 
direct production of pairs (p + 7 —> e^) and photo-pion- 
production (p + 7 —> ir). Figure 1 shows the energy loss rates 
P(p7 —> n/e^) of a proton in a black body radiation field by 
each of these processes for three different radiation temper- 
atures, computed from formulae given by Blumenthal (1970) 
and Stecker (1968). In the case of pair production, the approx- 
imation is made that most interactions occur near threshold, 
i.e., the energy of the incoming photon is neglected. A more ac- 
curate treatment has been given by Begelman et al (1990), but 
the differences (up to a factor of 2 for the highest temperature 
radiation field) are unimportant for our analysis. 

At any given temperature, there exists a value of the pro- 
ton energy above which losses by pion production are dominant 
and below which the losses are determined by pair production. 
The transition point between these two processes corresponds 
to a specific energy loss rate and, consequently, to a value of 
the electric field strength in the funnel. In stronger fields pro- 
tons lose energy mainly by pion production in weaker ones by 
pair production. The dependence of this critical value of the 
electric field E* , + in the funnel as a function of the radiation 
temperature is given by, 

Ei je± » 2.5 x 10" 8 • T 2 [Vcm" 1 ] (1) 

In Fig. 2 this relation is represented by the solid line. The 
Lorentz factor of a proton at which gains are balanced by losses 
determines the typical energy at which pairs are produced. In 
the region in which pion production dominates, pairs are pro- 
duced by the hard gamma-rays resulting from ir° -decay when 
they interact with a photon of the blackbody radiation field. 
The Lorentz factor of such pairs is about 100 times the pro- 
ton Lorentz factor and they decelerate not by inverse Compton 
scattering (ICS), but primarily by producing more pairs in the 
triplet pair production process (Mastichiadis 1991), thus form- 
ing an electromagnetic cascade in which few hard photons are 
present. This cascade continues until a bunch is formed (see 
below), since the ICS process is in the Klein-Nishina regime 
for funnel parameters appropriate to losses by pion production 
(see Fig. 2). One interesting consequence of the process of pion 
production is the conversion of the relativistic protons into rel- 
ativistic neutrons, which move ballistically and gain no further 
energy from the electric field. After only a few pion producing 
interactions, essentially all protons will have been converted 
into neutrons. If, as we assume, the radiation temperature in 
the funnel decreases outwards, the neutrons will not interact 
with thermal radiation in their path, but escape freely from 
the funnel, most of them decaying outside the central core of 
the AGN at a distance determined by their Lorentz factors 
(fl n [cm] « 3 x 10 13 - 7n ). 

If the interaction of protons occurs deep in the funnel where 
T and E { are below the line p7 —> ir/e^ then pairs are 
produced directly with Lorentz factors comparable to those 
of the primary protons. However, in contrast to the case of 
pion production, protons cannot convert to neutral particles 
and escape, but proceed to inject pairs throughout the whole 
funnel. The secondary pairs are decelerated in the funnel by 
Compton scattering until they achieve equilibrium either as 
single particles at 7 eq (for parameters below the line KN-T in 
Fig. 2) or as members of a bunch. 
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Fig. 1. The energy loss of protons P[eV/cm] per unit distance trav- 
elled in the funnel as a function of proton energy. (At saturation we 
have P = E* .) Shown are losses by pion production (dashed line) 
and e^ 1 pair production (solid line) in a black body radiation field 
with temperatures T = 10 5 " 5 K (l), T = 10 5 K (2) and T = 10 4 " 5 K 
(3). The horizontal dashed line defines the value of electric field E f 
(for T = 10 5 K) below which protons lose energy mainly by e 1 *- pair 
production and above which pion production dominates the losses 

2.2. Injection of electrons and positrons 

Let us assume that positrons of relatively low energy are in- 
jected into the electric field in the funnel. They are accel- 
erated by the component of the electric field parallel to the 
magnetic field until, at an energy 7 eq m e c 2 which depends on 
the local temperature of the radiation field T and the elec- 
tric field strength E l , their rate of energy gain is balanced by 
energy losses due to inverse Compton scattering. Synchrotron 
radiation, which can compete with inverse Compton losses for 
isotropically distributed particles if B>4.10 - T 2 is unimpor- 
tant for the accelerated particles, since these conserve the com- 
ponent of their momentum perpendicular to the magnetic field 
during acceleration. On scattering, the pitch angle is increased 
to at most a value ~ I/7. 

If the scattering takes place in the Thomson regime, i.e., 
if 7e q (3fcBT/m e c ) < 1, the equilibrium Lorentz factor (7 eq ) is 
given by 

7eq « 1.5 x 10 13 • (£ f ) 1/2 T- 2 . (2) 

In the conditions envisaged, positrons attain this Lorentz factor 
essentially instantaneously, after travelling a distance 

d eq « 10 19 (£ f )- 1/2 T- 2 cm , (3) 

which is much smaller than any macroscopic dimension of in- 
terest. 

Positrons coasting with Lorentz factor 7 eq in the Thomson 
regime transfer energy gained from the electric field directly 
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into high energy photons by ICS. However, in some parts of 
the funnel it is possible for the scattering to occur in the Klein- 
Nishina regime where 7 eq (3fcBT/m e c ) > 1. The line divid- 
ing these two parts of the parameter space (T,E ) is depicted 
in Fig. 2 and was obtained by setting the average energy of 
thermal photons in the reference frame of the electron (i.e., 
3fcBT7e q ) equal to the electron rest-mass energy. This gives 
the critical electric field -Ejcn-t as a mnc ti° n °f temperature: 



Ekn- 



1.6 x 10"°^ 
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Fig. 2. The separation of the different modes of interaction of protons 
and electrons in the narrow funnel of a thick accretion disk. The pa- 
rameter space is denned by the local values of the electric field in the 
funnel E* and the temperature of radiation X. The solid line py — Tr/e 
separates regions in which protons lose energy mainly by pion pro- 
duction (above this line) and e^ 1 pair production (below this line). 
The dashed line KN-T separates the regions in which single electrons 
move in the Klein- Nishina regime (above) and Thomson regime (be- 
low). Above the dash-dotted line E^{y / < 0.511 MeV) electrons 
comptonise thermal photons to energies > 0.51lMeV (these photons 
may produce additional e^ 1 pairs in 7 — 7 collisions). The two dotted 
lines give the range of parameters used in the models of CTA 102 
and QSO 1633+382 

It is also interesting to find out for which parameters (T,i£ f ) 
most of the photons produced in the comptonisation process 
have energies above 0.511 MeV. If such photons are created 
by both electrons and positrons, they are beamed in opposite 
directions, and are capable of producing additional pairs. 
We can estimate this condition easily by applying the formula 
for the average energy (E^) of scattered photons: 



(3fc B T)7e q , 



Replacing 7 eq in Eq. (5) by the value from Eq. (2), we find for 
the critical field E S11 (T) the relation 



£511 (T) « 6.2 x 10" 18 T 3 



(6) 



The corresponding critical line in parameter space is indicated 
in Fig. 2 by the dot-dashed line. 

An examination of Fig. 2 gives us an idea about the con- 
ditions in the funnel required for efficient production of 7-ray 
photons from inverse Compton scattering of low energy ther- 
mal photons by relativistic positrons. If electrons and positrons 
are injected at a location deep in the funnel where 7 eq lies in 
the Thomson regime, the 7-ray spectrum can be obtained sim- 
ply by integrating the contributions of positrons (the outgoing 
component) moving with Lorentz factor equal to the local value 
of 7e q throughout the funnel, provided the probability of inter- 
action between two 7-rays is sufficiently small. This restriction 
depends on both the luminosity of the source and the geometry 
of the funnel. For the case of blazars, we can make a simple esti- 
mate by taking a funnel of length r m ax = 10 cm and allowing 
for a flux of photons corresponding to an isotropic luminosity 
of L = 10 48 erg s -1 in gamma-rays of energy E y > 100 MeV. 
The photon number density in the funnel is then 



(7) 
(8) 



4irr max 2 cE 7 
2.10 8 cm" 3 



If we assume equal numbers of these photons travel inwards 
and outwards, the optical depth through the funnel to pair 
production in a head-on 7-7 collision is roughly 



< 2 x 10" 



(9) 
(10) 



where we have taken the interaction cross section appropriate 
to the interaction of two gamma-rays each of energy 100 MeV: 
<r 77 Ri 10 -28 cm 2 (Gould fc Schreder 1967). Thus, interactions 
between gamma-rays produced by inverse Compton scattering 
are unlikely in the blazar case. 

2.3. The Formation of Bunches 

In the case of injection of positrons in the Klein- Nishina regime 
(above the line KN-T), the situation is much more complicated, 
since 7-rays produced by Compton scattering can interact with 
soft photons from the funnel giving rise to an electromagnetic 
cascade. The sequence of events: 7-ray production by ICS —> 
pair production by the 7-ray —> acceleration of the pair to 
7e q —> 7-ray production by ICS, creates an avalanche of pairs 
similar to that discussed for the case of radio pulsars by Michel 
(1991). The Lorentz factor f CI corresponding to the boundary 
between the Thomson and Klein-Nishina regimes can be esti- 
mated as 



7cr 



m e c 
(3fc B T) 



2.10 a 



(11) 



(5) 



A soft photon of dimensionless energy a; s = hu s /mc 2 (with u s 
the frequency) which is scattered by an electron of Lorentz fac- 
tor 7 cr has its energy increased to roughly x = 47 2 r a; s /3 and 
is, therefore, just at the threshold for pair production with an 
unscattered photon, given by xx s > 2. Thus, any pair injected 
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in the Klein- Nishina region will be accelerated to an energy at 
which the inverse Compton gamma-rays it produces are capa- 
ble of producing more pairs. In contrast to the mean free path 
for pair production off another gamma-ray, the mean distance 
A 77 travelled before pair production with a soft photon is very 
short: 



(<T 77 n p h) 1 
8.10 23 



(12) 
(13) 



where n p h is the density of the soft photons and we have ap- 
proximated the interaction cross section by <r 77 « 3<tt/16, 
where <tt is the Thomson cross section (as is appropriate for 
interactions just above threshold). The secondary pairs are 
accelerated by the electric field, positrons moving alongside 
the seed particle and electrons in the opposite direction. Both 
types of particle attain a Lorentz factor 7 eq within a distance 
deq given by Eq. (3) and subsequently produce photons capable 
of creating the next generation of pairs. This cascade continues 
as long as the particles move in the unmodified external electric 
field. However, once the accumulation of positive charge around 
the initial positron is sufficient to screen the external field from 
the particles immediately behind it, negatively charged parti- 
cles created in the bunch will no longer be removed. Although 
the process of pair production by each electron and positron 
in the bunch continues, increasing the total number of parti- 
cles present, and thus the total drag force exerted by the soft 
photons, the net charge of the bunch ceases to rise, reducing 
the equilibrium value of the Lorentz factor at which the bunch 
moves. Eventually, the bunch reaches a critical value of the 
Lorentz factor below which pair creation is no longer possi- 
ble. Given constant external parameters, we expect the bunch 
to remain close to this critical Lorentz factor: it cannot move 
much faster, since then more pairs are created increasing the 
drag, and it cannot move much slower, since the number of 
particles would then remain constant, and with it the accel- 
eration and drag forces. This scenario closely resembles that 
described by Sturrock (1971) and Michel (1991) for the for- 
mation of dense almost neutral bunches of pair plasma in a 
pulsar magnetosphere. However, in the case of pulsars, the role 
of the background thermal radiation is taken over by virtual 
photons associated with the strong magnetic field. 

In order to estimate the total number of pairs in a bunch 
we follow Michel (1991) and consider a single positron which 
crosses the KN-T line in the parameter space of Fig. 2. At each 
step of the cascade, particles are produced which move almost 
precisely along the electric field direction at almost exactly the 
speed of light. The spatial dispersion of particles along the field 
direction is, therefore, extremely small, and the bunch develops 
into a pancake or disk whose extension perpendicular to the 
field direction is caused by the small angle (~ l/~f C i) which 
scattered photons subtend with the direction of the electric 
field. For each fresh generation, therefore, the new electric field 
line on which pairs are born is displaced from the old one in a 
random direction (perpendicular to the field) by a distance 



l/(<TTre p h7cr). . 



(14) 



After N cascade generations, the random walk across electric 
field lines produces a disk of radius approximately AJV 1 '' 2 , and 
area S given approximately by 



The minimum number of positrons ra+ required to form a 
bunch can be estimated assuming all electrons are removed 
until the field produced by the positron disk exactly balances 
the externally applied field just behind the disk. Using Gauss' 
law we find: 



SE f /(4Tre) 
10 34 E'N/T 1 . 



(16) 



The development of the cascade is described in Fig. 3. The 
number of positrons ra+ (N) in the leading front after N inter- 
actions is given by the recursion relation 



n + (N) = n+(N - 1) + n + (JV - 2), 



(17) 



together with the initial values ra+(0) = 1, ra+(l) = 1. For large 
N we have approximately 



n + (N) « exp[(JV-3)/2]. 



(18) 



For a reasonable electric field strength E l = 4 x 10 3 V cm -1 and 
temperature T = 5 x 10 5 K we find, by comparing equations 
(16) and (18), that saturation sets in at about N « 80 at which 
point the front contains some 5 x 10 16 positrons (see Fig. 4). 



Avalanche Development 



Time 




1.8 X 10 28 JV/T 4 



(15) 



Fig. 3. A simplified one-dimensional picture of avalanche develop- 
ment. An interaction (Compton scattering or 7—7 pair production) 
is assumed to occur after a particle (e 1 * 1 or photon) has travelled ex- 
actly one interaction length. The electric field maintains the speed 
of charged particles at approximately that of light, with electrons 
moving to the left and positrons to the right 

At this point, our simple one dimensional picture has two 
fronts, one containing 5 x 10 positrons moving outwards and 
another containing an almost equal number of electrons moving 
inwards. Between them, there is a fairly smooth distribution of 
charged particles which, however, find themselves in a region of 
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Fig. 4. The almost exponential growth of the number of positrons in a 
bunch as a function of the number of interactions N (solid line). The 
dashed line depicts the number of positrons required in the bunch 
in order to cancel the external electric field after N interactions. 
The point of intersection gives the estimated bunch dimension. The 
parameters used are 4 X 10 V cm" _1 and T = 5 X 10 5 K. 



zero electric field. In reality, however, the finite transverse ex- 
tent of the disk remains very small compared to an interaction 
length, as can be seen from Eq. (15), so that the shielding effect 
of the leading disk will not extend to following disks, but will 
merely ensure that negatively charged particles subsequently 
created close to the leading front are no longer removed from 
it. 

Models very similar to this one involving bunches of 
charged particles have been discussed by several authors in 
connection with the radiation mechanism of radio pulsars (for 
a review, see Melrose 1992). There are two main problems with 
them: (i) whether or not the bunches are stable and (ii) whether 
or not they are capable of emitting radio waves coherently, as 
required to explain the radio pulses. Only the first of these 
points is relevant here. The stability of bunches is, of course 
closely related to the process responsible for their formation. 
Early work suggested that the growth of waves in the two- 
stream or beam-plasma instability could lead to a series of 
charged bunches. However, these are rather quickly destroyed 
in the nonlinear saturation phase of the instability. Only in the 
case of a beam penetrating a completely charge separated back- 
ground (Elsasser fc Kirk 1976) does it seem possible to have 
bunches which exist for many plasma periods (i.e., many times 
the inverse of the plasma frequency). Even so, there are diffi- 
culties connected with the coherent radiation of such bunches 
(Asseo et al 1980). 

Quite a different approach is required to the problem of 
bunch stability in the model proposed by Michel (1991). Here 
it is assumed that regions around a pulsar exist in which a 
strong electric field oriented parallel to the magnetic field main- 
tains a vacuum gap. In close analogy with the case we describe 
above, bunches are created out of the vacuum by interaction of 
a photon with the magnetic field, and the process of curvature 
radiation by the created pair provides more photons to con- 
tinue the cascade. The particle energy is limited by radiation 
reaction, and the number of particles in a bunch exponentiates 



in time until the oppositely charge lepton starts to be bound to 
the bunch. Plasma instabilities are not thought to play a major 
role in the evolution of the cascade, but, instead, work has con- 
centrated on the dynamics of the individual charged particles 
under the influence of radiation drag and the ambient electric 
field. Currently, there are two points of view: Fawley (1978), 
proposes that an electron positron pair injected into such a 
gap produces two 'flame-fronts' moving in opposite directions 
along the magnetic field, rather like the outward and inward 
propagating electron and positron bunches described above. 
However, on the basis of one- dimensional simulations, he hy- 
pothesises that the region between the flame fronts has zero 
electric field and is filled with an almost neutral plasma. On 
the other hand, Michel (1991) proposes a picture in which also 
the interior region between the flame fronts organises itself into 
bunches. It is this picture which we adopt. The physical jus- 
tification lies in the inherently two-dimensional nature of the 
problem: the transverse dimension of a flame front can be esti- 
mated (cf. Eq. 15) to be very much smaller that the separation 
of the fronts. Thus, the trailing cloud of electrons and positrons 
which are completely shielded from the ambient electric field 
in the one- dimensional picture must be almost completely ex- 
posed to it in reality. It is thus plausible that they also evolve 
in such a way as to create their own flame fronts, i.e., bunches. 

2.4- The Inverse Compton Emission of Bunches 

Although the stability of the proposed bunches remains to be 
demonstrated, it is interesting to analyze the fate of a bunch 
once it has accumulated sufficient charge to bind its electrons. 
If we define fb{n+ ,"-) to be the equilibrium Lorentz factor 
of a bunch containing n_ electrons and ra+ positrons, we can 
balance gains and losses to find: 

rjeE 1 = -Pics(7b,T), (19) 
c 

where 

r, = (n+ -n_)/(n+ +n_), (20) 

Pics{"fb,T) is the energy loss rate a single electron or positron 
with Lorentz factor 7b by inverse Compton scattering of black 
body photons with temperature T, and — e is the charge of an 
electron. When saturation occurs, the value of the parameter 
T) should not be too different from unity: j/o = J?(r m i n ) « 1. 
However, the process of pair production continues to cause a 
rapid rise in the number of pairs and consequent decrease of j/, 
so that we expect the bunch to settle to an equilibrium Lorentz 
factor just inside the Thomson regime: 7b ~ 7 cr - We can then 
use Eq. (19) to determine j/, given that the net charge of the 
bunch is roughly the number of positrons it contained when 
electrons started to be bound to it. 

The subsequent evolution of the bunch through the funnel 
is determined by the electric field and temperature profiles. 
Two cases can be distinguished: 

1. If E l T~ 2 increases along the path, ij decreases steadily and 
is given by Eq. (19). The Lorentz factor of the bunch is then 
always equal to f CI . 

2. If, on the other hand, E T~ 2 decreases, the Compton drag 
decelerates the bunch below f CI , pair creation stops, and r] 
remains constant. 
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In each case, we can estimate the spectrum radiated by a 
bunch using the 5-function approximation to Compton scatter- 
ing, in which the number F of photons of frequency v emitted 
by a single particle per second per unit frequency interval is 
given by: 



dF 
dvdt 



4<T T C7 2 i7 ra d 



8(v - fo), 



(21) 



where f7 ra d is the energy density in the background photons and 
uo = 4(3A;bT)7 2 /(3/i) is the frequency of scattered radiation. 
Once created, a gamma-ray photon leaves a bunch freely, since 
the optical depth to scattering through the bunch, ji+ctt/ S, is 
very small (typically ~ 10" 15 ), as can be seen from Eqs. (15) 
and (18). The total emission per unit solid angle dJ- / di/dQdt 
is found by integrating over the path of a bunch: 



dJF 



di/dfidt 



dr JV t 



dF 
' dvdt 



(22) 



where N is the rate of injection of bunches per unit solid angle, 
JVb = iVoJjo /v is the total number of electrons and positrons in 
the bunch, rjoNo = (ra+ — n_) is the net bunch charge (which 
remains constant after bunch formation) and No = Nb at the 
moment of bunch formation. The emission of a single particle is 
beamed into a very small angle (w l/7eq) about the direction 
of motion. However, even a slight flaring of the funnel walls 
should cause the electric field lines and hence the direction of 
motion of the particles to diverge, allowing an observer to sam- 
ple the entire cone of single particle emission. Consequently, it 
is appropriate to use the expression for the emissivity inte- 
grated over all emission angles in Eq. (22). Adopting power 
law profiles according to: 



V '"min / 

T = T max (^-)~ q 

\ * min / 

we find the spectra 



(23) 



'2.4 x 10 14 (3.3 x lO 29 )^- 1 ^^- 1 

XjVJVoJJormax-Bmax 



dF 



di/dfidt 



T -(j>-l)/« -2-(p-l)/ ? 

in case 1 (p < 2g) 



(24) 



2.4 x 10 14 (1.9 x lO 37 )^- 1 ^ 39 -^ -pi" 1 

vT -3(p-l)/(3 ? -p) -2-(p-l)/(3 ? -p) 

in case 2 (p > 2g) 



It is interesting to note that flat spectra (of index < 2) 
can be found for p < 1 in case 1, provided the temperature 
decreases outwards (g > 0). In case 2 flat spectra are also 
obtained for p < 1, provided the temperature gradient is suffi- 
ciently steep (g > p/3). 

3. Application to 7-ray emission from blazars 

In order to investigate the applicability of a bunch scenario 
to 7-ray emitting AGNs, we have fitted the 7-ray spectrum of 
two blazars (CTA 102 and QSO 1633+382), recently observed 
by EGRET (GRO) at energies above 100 MeV. CTA 102 is a 



7-ray emitting AGN with one of the softest spectra so far ob- 
served by EGRET (spectral index 2.6 ±0.2, Nolan et al. 1993). 
In contrast, QSO 1633+382 shows one of the hardest spectra 
in the EGRET sample (spectral index 1.86 ± 0.07, Mattox et 
al. 1993). The strategy we adopt to obtain a fit is first to fix the 
maximum and minimum radii and temperatures, according to 
ideas concerning the structure of thick accretion disks and to 
observational indications of the likely temperature range. The 
parameter q is thereby determined (see Table 1). We then have 
the two choices indicated in Eq. (24) for obtaining a value of 
p which gives a good fit to the spectrum. In case 2, in which 
a bunch moves with Lorentz factor less than f CI (defined in 
Eq. 11) CTA 102 would require p = 1.8, which is in contra- 
diction with the condition p > 2g, needed for bunch emission. 
However, this case applies not only to bunches of fixed charge, 
but also to individual particles, provided they scatter in the 
Thomson regime. In fact, the required spectrum can be ob- 
tained if the particles remain beneath the KN-T line in Fig. 2 
and do not form bunches i.e., if the electric field within the fun- 
nel is sufficiently small (£ miI < 7x 10 V cm - for the assumed 
range of temperature). If, on the other hand, the parameters 
of the funnel imply that particles cross the KN-T, then they 
immediately form bunches, and, since p < 2g, we must employ 
the formula appropriate to case 1. Since a bunch contains a 
large number of particles, that part of the emission produced 
above the KN-T dominates any contribution from the single 
particle whilst it moved below the line. The value of p required 
for case 1 is 1.5, consistent with the requirement p < 2g. 

The situation for the hard spectrum source QSO 1633+382 
is similar. Once again, a straightforward application of the for- 
mula pertaining to case 2 results in p = 0.19, in conflict with 
the condition p > 2g. A fit is possible for individual particles 
scattering in the Thomson regime (i.e., below the KN-T line in 
Fig. 2) only if £ mix < IV cm -1 . Consequently, we concentrate 
here on case 1. 
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Fig. 5. The 7-ray spectrum of the blazar CTA 102 (z % 1.) observed 
by EGRET (Nolan et al. 1993) and the fit obtained in the bunch 
scenario, with the parameters given in Table 1. 



8 



% -8 



\ -10 

C 
o 
+-> 

o 

X 
D. 

M -12 
o 



-14 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
QSO 1633+382 


: K 

\ 1 


- 


\ 




\ 

\ 

\ 

— \ 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


\ 

\ 

\ 

\ 

\ 

\ 



Table 1. Model parameters 



2 3 4 

log E 7 [MeV] 



Tig. 6. As in Fig. 5, but for the blazar QSO 1633+382 (z « 1.8) 
(Mattox et al. 1993) 



In calculating the spectrum we relax the delta function 
approximation to Compton scattering in favour of the for- 
mula given by Blumenthal and Gould (1970; Eq. 2.48). The 
fits to the 7-ray spectra of CTA 102 and QSO 1633+382 are 
presented in Figs. 5 and 6. They are obtained with the tem- 
perature and electric field profiles given by Eq. (23), using the 
parameters shown in Table 1. These profiles are also marked in 
Fig. 2; they correspond to emission by bunches of charge. The 
fits obtained for emission by single particles in the Thomson 
regime are not shown. They differ only in that the power-law 
extends to lower frequencies without the flattening visible in 
Figs. 5 and 6. 

For the models listed in Table 1, we estimate, using 
Eqs. (16) and (18) and the values of T max and the num- 

ber of positrons initially in the bunch to be roughly 10 15 for 
each source. Then, normalising the model flux to that observed 
by EGRET, we use Eq. (22) to find the rate of injection of 
bunches: 



N 



J3.10 22 
|3.10 21 



bunches sr 
bunches sr~ 



for QS01633 + 382 
for CTA102 



(25) 



where we have assumed the cosmological parameters f2 = 0.5 
and Ho = 75km/s/Mpc. Assuming an opening angle of 5° for 
the funnel, Eq. (25) implies a t otal of 2.10 20 and 2.10 19 bunches 
per second for QS01633+382 and CTA102, respectively. 

Figure 7, shows the evolution of the Lorentz factor and 
particle content of a bunch in the profiles used for modelling 
the two blazars. For each source, the evolution follows case 1, 
resulting in a steady increase in the Lorentz factor outwards. 
The number of particles in the bunch JVb increases too, keeping 
the scattering just on the boundary between the Klein- Nishina 
and Thomson regimes. 

It is important to note that the model we discuss is able to 
fit reasonably the broad range of spectral indices in 7-ray en- 
ergy range reported from AGNs by EGRET. In general, harder 
spectra are obtained when a large part of the funnel has a 
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Fig. 7. (a) The Lorentz factor 71, and (b) the number of electrons and 
positrons in the bunch as a function of position in the funnel for 
the parameters of CTA 102 (solid line) and QSO 1633+382 (dashed 
line) 



relatively strong electric field and low temperature. The low 
frequency inverse Compton emission arises in the inner parts 
of the funnel. Since the emitting particles are accelerated es- 
sentially instantaneously to 7 eq , and because they leave the 
funnel and its field of soft photons without being allowed to 
cool, there is a lower cut-off in the energy spectrum of radiating 
particles, which leads to a turn-over in the photon spectrum 
at about 10 MeV. The inverse Compton spectrum at lower fre- 
quencies eventually becomes completely flat (Blumenthal fc 
Gould 1970). However, we do not attempt to explain the spec- 
trum in the hard X-ray region within our model. All the 7- 
ray spectra in our model show a cut-off at few tens of GeV. 
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This is determined by particles which emit close to the Klein- 
Nishina/Thomson transition in the outer parts of the funnel, 
where the temperature is of the order of 10 K. In one case, 
Mkn 421, 7-ray emission has been reported at TeV energies 
(Punch et al. 1992). This source is, however, not typical of the 
EGRET sample (it is one of the closest objects detected so 
far by EGRET and has a significantly lower luminosity than 
observed from other AGNs), and may operate by a different 
radiation mechanism. In a model such as ours, TeV photons 
do not arise naturally, since this would require either a much 
lower temperature in the radiation field or a strong anisotropy 
(Burns fc Lovelace 1982). One possibility - that injected pro- 
tons are converted into neutrons which decay at large distance 
from the core - can be discounted: this requires proton injection 
above the solid line in Fig. 2 in order to convert the charged 
particle into a neutron. However, the associated positron from 
pion and muon decay is injected in a region above the KN-T 
line, and so will produce a bunch containing a large number of 
positrons (see Eq. (16)). The energy extracted from the elec- 
tric field by these positrons would by far exceed that taken off 
by the neutron. 

4. Summary and conclusions 

We investigate the possibility that the 7-ray emission recently 
detected from more than 20 AGNs is produced in the nar- 
row funnel of a thick accretion disk by relativistic electrons 
and positrons which comptonise soft thermal photons from the 
funnel walls. These particles, injected directly or as secondaries 
from hadronic interactions, are assumed to be accelerated in a 
longitudinal electric field which may be induced in the inner 
parts of the funnel along magnetic field lines which thread the 
horizon of a rotating supermassive black hole. Other scenar- 
ios which produce aligned magnetic and electric electric fields 
(such as reconnection in an MHD jet) produce basically similar 
effects, but the details of the spectrum may be different. As a 
result of the balance between energy losses by inverse Comp- 
ton scattering and energy gains from electric field, electrons 
and positrons move through the funnel in either the Thom- 
son or Klein-Nishina regimes of scattering (see Fig. 2 for the 
corresponding parameter space in T and E l coordinates). In 
the Thomson regime, the electrons move with the equilibrium 
Lorentz factor defined in Eq. (2) and efficiently comptonise soft 
photons into the 7-ray energy range. 

The propagation of electrons in the Klein-Nishina regime 
is much more complicated. We argue that single electrons pro- 
duce dense bunches of plasma containing relatively low net 
electric charge. A bunch, on its subsequent path through the 
funnel, may grow by producing captive pairs or decelerate into 
the Thomson regime, according to the T and E profiles in the 
funnel. Our bunch production scenario resembles the produc- 
tion of bunches of plasma in pulsar magnetospheres pro- 
posed by Michel (1991). 

A bunch moving through the funnel extracts energy from 
the electric field and converts it via inverse Compton scatter- 
ing into 7-ray photons. The Lorentz factor of the bunch is 
self-regulated such that it is either close to or below the value 
at which the scattered photons could produce pairs on inter- 
acting with the thermal photons. As a result, softer 7-rays are 
produced in regions where T is higher (in the inner part of 
the funnel) and harder 7-rays in regions of lower T (the outer 
parts of the funnel). We model the 7-ray emission from blazars 



in a bunch scenario (sect. 3), fitting, as examples, the 7-ray 
spectra of two blazars (CTA 102, Fig. 5 and QSO 1633+382, 
Fig. 6). However, it is also possible to obtain fits to the spectra 
without forming bunches. This requires radiation by single par- 
ticles which scatter the soft photons in the Thomson regime. In 
order that pair creation be avoided, a relatively weak electric 
field is necessary (see section 3) - for the hard spectrum source 
QSO 1633+382 this condition is particularly restrictive. In our 
model, the spectrum below the 7-ray region flattens, so that 
we do not produce the nonthermal X-ray and optical emission. 
If we were to increase the radiation temperature in the inner 
parts of the funnel, X-ray photons would result, but this would 
imply a thermal component of much higher temperature than 
usually observed in the soft X-ray excesses of AGNs. Instead, 
we propose that the nonthermal X-rays and optical photons 
are emitted outside the funnel either by comptonisation of the 
anisotropic radiation from the funnel, or by synchrotron radia- 
tion. A basic problem encountered by all models of gamma-ray 
emission from blazars is that of absorption by the hard X-ray 
component. Since our source of gamma-rays is close to the 
central object, an isotropic distribution of hard X-ray photons 
would absorb it completely. This problem can be avoided only 
by assuming the nonthermal X-rays are beamed in the same 
direction as the gamma-rays. A beaming angle of 5°, for ex- 
ample, implies that a 100 MeV gamma-ray requires a target 
photon of 1 MeV in order to pair produce, whereas a 5 keV 
target suffices in the isotropic case. From Eqs. (7) and (9), one 
then finds that the effect of absorption is negligible. 

If protons are injected deep in the funnel into the accel- 
eration region, their interaction with soft photons is by direct 
production of pairs and/or it's (which then decay, produc- 
ing pairs). The relative importance of these two processes 
depends on the T and E in the injection region (these two 
regions are separated in Fig. 2 by the line pf —> ir/e^). There 
are significant differences between direct injection of protons 
and electrons in the acceleration region: 

1. The pairs, produced as a secondaries by the decay of it's, 
have such high Lorentz factors that the process of triplet 
pair production determines the bunch formation. However, 
the fate of the bunch is similar to the case of direct injection 
of low energy electrons. 

2. The interaction of protons in the regime of direct pair 
production (below the line pf —> ir/e^) produces an injec- 
tion of pairs along the whole path of proton's motion in the 
funnel. If this happens below the line -B 7 (> / < 0.5MeV), 
the pairs are slowed down to the Thomson regime and 
comptonise the soft thermal photons into the hard X-ray 
energy range. If the injection of pairs occurs above this line, 
the 7-ray photons, produced in the comptonisation process, 
can, in principle, interact with each other to produce pairs 
(electrons and positrons produce beams of 7-rays moving 
in opposite directions). The timescale of this instability 
depends on the optical depth in the funnel. If it is short 
enough, the entire funnel is filled with dense low energy 
plasma which may partially annihilate (producing broad 
annihilation feature) and may be partially expelled outside 
the funnel by radiation pressure (a mildly relativistic 
jet may be created). In blazars, this optical depth is very 
small, but the scenario could conceivably be applied to the 
recently reported transient source of annihilation radiation 
in the Galactic Center 1E1740. 7-2942 (Bouchet et al. 1991; 
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Sunyaev et al. 1991). However, we do not pursue this line 
of investigation here. 
3. Almost all of the protons, which interact with soft pho- 
tons in the region characterised by parameters T and E 
above the line p7 —> ir/e^ , are converted after a few inter- 
actions into relativistic neutrons. The neutrons can leave 
the funnel without significant energy losses because they 
are below the threshold for ir production in the outer re- 
gions of the funnel (T and E decreases towards the outer 
parts of the funnel). They decay at a distance determined 
by their Lorentz factor. However, a much larger luminosity 
is expected of the bunch produced by the positron resulting 
from decay of the associated ir + . 

Since plasma escaping from the funnel carries a net elec- 
tric charge, it may be supposed that the central black hole 
(and/or inner parts of the accretion disk) will accumulate an 
opposite charge after certain time, preventing the stationary 
operation of the picture discussed here. However, as noted by 
Lovelace (1976) this charge can be neutralised by a flow of 
current through the accretion disk. 

As showed in Figs. 5 and 6, our model can provide high 
energy 7-ray spectra with a the broad range of spectral in- 
dices. However, the spectra show a characteristic sharp cut-off 
at an energy of a few tens of GeV. This is determined by the 
temperature in the outer parts of the funnel. It is difficult to 
determine at present, on the basis of CGRO observations, if 
such a cut-off is indeed present in the spectra of AGNs. One 
clear counter-example is Mrk 421 which has been detected at 
TeV energies. However, this source is not a typical member 
of the sample of powerful 7-ray emitting AGNs reported by 
EGRET. 
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